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Bound-state eigenenergy outside and inside the continuum 
for unstable multilevel systems 
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The eigenvalue problem for the dressed bound-state of unstable multilevel systems is examined 
both outside and inside the continuum, based on the level Friedrichs model which describes the 
couplings between the discrete levels and the continuous spectrum. It is shown that a bound- 
state eigenenergy always exists below each of the discrete levels that lie outside the continuum. 
Furthermore, by strengthening the couplings gradually, the eigenenergy corresponding to each of 
the discrete levels inside the continuum finally emerges. On the other hand, the absence of the 
eigenenergy inside the continuum is proved in weak but finite coupling regimes, provided that each 
of the form factors that determine the transition between some definite level and the continuum 
does not vanish at that energy level. An application to the spontaneous emission process for the 
hydrogen atom interacting with the electromagnetic field is demonstrated. 



PACS numbers: 42.50.Vk, 42.50.Md 



I. INTRODUCTION 

The theoretical description of the unstable quantum 
system often refers to the system of a finite level coupled 
with the spectral continuum. In weak coupling regimes, 
the initial state localized at the finite level undergoes 
exponential decay pj. However, by changing the cou- 
plings to stronger regimes, instead of the total decay a 
partial one can occur Q. This means that the superpo- 
sition between the states localized at the finite level and 
the continuum forms the dressed bound state, that is, a 
bound eigenstate extended over the total Hilbert space. 
The formation of the bound eigenstate is of great inter- 
est in the study of the various systems having to do with 
such matters as the photo detachment of electrons from 
negative ions @, Q and the spontaneous emission of pho- 
tons from atoms in photonic crystals 0, IE • It is then 
clarified that the energy of the bound eigenstate depends 
not only on the strength of the couplings but also on the 
relative location between the electron bound-energy and 
the detachment threshold 0, 01 > or between the energy 
of the atomic frequency and the continuum edge of the 
radiation frequency [|J, Q ■ Further research has been di- 
rected to those eigenstates, aiming at the decoherence 
control HEI 

In these analyses, however, single-level systems are 
treated often, while multilevel systems are less exam- 
ined. In the latter, some peculiar time evolutions are 
theoretically observed: steplike decay , decaying oscil- 
lati on jrij . and various long-time nonexponential decays 
[T3. Il3|. These peculiarities are never found in single- 
level approaches. Furthermore, to the author's knowl- 
edge, the possibility of a bound-state eigenenergy "in- 
side" the continuum has not been studied except in a 
special multilevel case where all form factors are assumed 



to be identical |1[T^|. 

In the present paper, we attempt to examine the eigen- 
value problem for the dressed bound-state in multilevel 
cases, based on the AT-level Friedrichs model IT^| . 
allowing some class of form factors, including identical 
cases. We show that for the discrete energy levels ly- 
ing outside the continuum, the bound-state eigenenergy 
always remains below each of them. Moreover, by in- 
creasing the couplings, the bound-state eigenenergy cor- 
responding to each of the discrete levels inside the con- 
tinuum can emerge out of that continuum. For the 
bound-state eigenenergy inside the continuum, we can 
only prove its absence in weak coupling cases under the 
condition that the form factors do not vanish at the en- 
ergy of each level. This result is just an extension of the 
lemma already proved for a system with identical form 
factors 14]. An upper bound of the coupling constant for 
the case of no bound-state eigenenergy being inside the 
continuum is also obtained explicitly. We apply this re- 
sult to the spontaneous emission process for the hydrogen 
atom under the four-level approximation. 

In the next section, we introduce the TV-level Friedrichs 
model and its eigenvalue problem. In Sec. IIIII we con- 
sider the eigenvalues outside the continuum, with resort 
to the perturbation theory about the eigenvalue of her- 
mitian matrix. The discussion developed here helps us 
to undertake the problem for the inside case, which is ar- 
gued in Sec. IIVI Concluding remarks are given in Sec. [V] 
We also present an appendix where both the small and 
large energy behaviors of the energy shift are studied in 
detail. 



II. THE A-LEVEL FRIEDRICHS MODEL AND 
THE EIGENVALUE PROBLEM 
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The iV-level Friedrichs model describes the AT-level 
system coupled with the continuum system. The total 



2 



Hamiltonian H is defined by 



H = H + XV, 



(1) 



where A € R is the coupling constant. We here define the 
free Hamiltonian Hq as 

N 

Hp = y~lw n \n)(n\ + / cj|cj)(w|p(Lj)da;, (2) 

„_i 



By setting this into Eq. (|3J), we have the normalization 
condition 



o 



(i 



u>-E- 



du> < oo, (9) 



which is the essential of the localization of dressed bound- 
state. 



where it was assumed that wi < loi < ■ • • < Wjv. \n) and 
\ui) satisfy the orthonormality condition: (n\n') = 5 nn >, 
(u)\u>') — S(lo — lo')/ p(uj), and (n\aj) = 0, where 5 nn > is 
Kronecker's delta and 5(ui — lo') is Dirac's delta function. 
p(lu) is a nonnegative function interpreted as, e.g., an 
electromagnetic density of mode, and ft = {lo\p(lo) ^ 0} 
is a specific region, like the energy band allowed by the 
electromagnetic mode. The interaction Hamiltonian V 
describing the couplings between \n) and \lo) is 

N 

V = J2 [v n (u)\u)(n\+v*{u)\n){u\]p(u)du, (3) 



where (*) denotes the complex conjugate and v n (u>) is 
the form factor characterizing the transition between |n) 
and \u>). We here assumed that v n S L 2 (0, oo), i.e., 



\v n (cu)\ 2 p(uj)duj < oo. 



(4) 



For clarity of discussion below, we assume that p(lo) = 1 
for u> > and otherwise, so that = [0, oo). Then, we 
merely write Jq by / °°, and the outside of the contin- 
uum means the half line (— oo, 0). An extension of fl to 
more general cases, such as gap structures, is not diffi- 
cult, however, our facilitation could extract the essential 
of the matter. 

Let us next set up the eigenvalue problem for this 
model. We suppose that the eigenstate corresponding 



to the eigenvalue E is of the form \ue) = E«— 1 
J f '(lo)\lo) du> , and it is normalizable, i.e., |17| 



Cn\n) + 



\c n \ 2 + \f(iv)\ 2 div< X> . (5) 



N 

(u E \u E ) = Ic„' 2 

n=l 



Then, the eigenequation H\ue) — E\ue) is equivalent to 
the following ones, 

/>oo 

io n c n + X v*(u)f(uj)dw = Ecn, Vn = 1,...,N, (6) 
Jo 



N 



J f(cj) + X^2c n v n (tj) = Ef(iv) 



n=l 



Equation J7J) immediately implies 



to-E 



/H = -a 



(7) 



(8) 



III. BOUND-STATE EIGENENERGY OUTSIDE 
THE CONTINUUM 

We first review the results on the negative-eigenvalue 
problem for N — 1, the single-level case [l^. If E < 
0, the integral in Eq. I© always converges under the 
condition In fact, 



K-il 



i«.Mi\ d „<ari,,Mi^<», 



l^ + |si| 2 



E 



(10) 

Thus, the substitution of f(u>) into Eq. I© is allowed. 
By introducing the function n(E) as 



k(E) = ui - X 2 



■du), 



Eq. © reads [jj 



k(E) = E, 



(11) 



(12) 



which is either an algebraic or transcendental equation 
of E, depending on v\(lo). k(E) has two important prop- 
erties as follows, 



k(E') > k(E), and k(E) < to u 



(13) 



for all E and E' satisfying E 1 < E < 0. The former 
means that k(E) is monotone decreasing in E. Therefore, 
there is only one solution (negative eigenvalue) E of Eq. 
(fT2*|) if and only if 



lim k(E) = u)\ — lim A 2 / 



LO — E 



■du < 0. (14) 



When E > 0, E should be a zero of Ui(cj) so that Eq. 
© holds. This is discussed in detail in Sec. IIVI 

Let us now turn to the TV-level case. Corresponding to 
Eq. (|10J) . this time we have that 



Er^lCr^nHI 2 ^ < En=l / °° KH| 2 G^ ^_ ^ 



lw+ \E\ 



\E\ 



(15) 

and Eq. J3J is satisfied again, where we used that 
E^=i l c "l 2 ^ 1- Substituting Eq. © into ©, one ob- 
tains 

N 

^ [w„<W - A 2 ^'^)]^' = E 1 ^, (16) 
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where 



Snn' { z ) 



v*(aj)v n > (uj) 



du>, 



(17) 



with z S C\[0, oo). For a later convenience, we here 
introduce an N x N matrix S(z) with the components 
S nn >(z). Note that S(E) for E < turns out to be a 
Gram matrix |20j , which is positive semidefmite. One 
obtains the following property of S(E): 

Lemma III.l S(E') < S{E) for E' < E < 0. 

Proof : We have that 

S nn >(E) - S n ,AE') = {E-E')T nn ,{E,E'), (18) 

for all E and E' satisfying E' < E < 0. We here intro- 
duce the matrix T{E, E') whose components are 



T nn '(E, E) := / 
Jo 



v*(uj)v n ,{iv) 
(ui — E){u) — E') 



duj. 



(19) 



Note that since T(E, E') is a Gram matrix, it is positive 
semidefmite. Therefore the proof is completed. D 

We also introduce the matrices Kq and K(E) — Kq — 
X 2 S(E) with components 



and 



K nn '{E) :— u! n 5 nn i — X 2 S nn '(E), 



(20) 



(21) 

respectively. For any E < 0, K(E) becomes a hermitian 
matrix, and thus there are A eigenvalues of K(E). We 
denote them by {K n {E)}^ =l , where n\{E) < k 2 (-E) < 
. . . < kn(E). The existence of a nontrivial solution {c n } 
of Eq. I |lt>[l is guaranteed if and only if there exists a 
negative E to satisfy 

K n (E)=E, (22) 

for a certain integer n. As in the former part of Eq. (|13fl . 
K n (E) has the following property: 

Lemma III. 2 For any fixed n, K n {E') > n n (E) for 
E' < E < 0. 

Proof : We see from Eq. (|T%)l that 

K(E) ~ K(E') =-(E- E')\ 2 T(E, E') < 0, (23) 

for E' <E < 0. Then, by using the Theorem 4.3.1 in Ref. 
[20j . the following inequality between the eigenvalues of 
K(E), K(E'), and T(E, E') holds [HJ, 

K n (E') - (E - E')X 2 t n (E,E') 

< K n (E) < K n (E') -{E- E')X 2 t 1 (E,E'), (24) 

where T n (E, E') denotes the n-th eigenvalue of T(E, E'). 
Note that since T(E, E') > 0, all T n (E 7 E') > 0. Then, 
— (E - E')ti(E, E') < for E > E', and the inequality 

K n (E) < K n (E'), (25) 

immediately follows from the last part of Eq. I|24|) . D 

We also have the statement below, which corresponds 
to the latter part of Eq. (|13fl . 



Lemma III. 3 For any fixed n, K n (E) < uj n for all 
E < 0, and lim n n {E) — uj n . 



_E->-c 



Proof: From Eq. (gTJl and Theorem 4.3.1 in Ref. |2C 
again, one obtains that 



u n - X 2 a N (E) < K n (E) <0J n - X 2 ai{E), 



(26) 



where cr n (E) denotes the n-th eigenvalue of S{E). If we 
recall the fact that S(E) > implies cr n (E) > for every 
n, the above inequality reads 

< X 2 a 1 {E) <uj n - K n (E) < X 2 a N {E). (27) 

Asymptotic behavior of the right-hand side of the above 
can be evaluated from Eq. l(T7|) as 



N 

( £ )<«k S(£ ,,< m x:/ 



\v n {u)\ 2 dcu -> 0, (28) 



as E — > — oo, and thus the lemma is proved. D 

Therefore, summarizing Lemmas 1111.21 and MI. 31 we 
obtain 

Theorem III. 4 If liniETO K n (E) < up to n — M, 
then each of the K n (E) for n = 1, . . . , M intersects E 
only once, so that M negative eigenenergies of H exist. 
In particular, if Hq has AL. negative eigenenergies, i.e., 
u> n < up to n = A_ , then N_ negative eigenenergies of 
H, denoted by E n , exist and satisfy E n < oj n . 



We also see from Eq. (|27?)) that 

(E)<uj n -X 2 a 1 {E). 



(29) 



This means that when |A| is large enough, every K n (E), 
even originating from a positive w n , becomes negative, 
unless o~\{E) — 0, i.e., the v„(w)'s are linearly dependent 
[20|. More precisely, the following statement holds. 

Proposition III. 5 Suppose that only N- ln d form fac- 
tors are linearly independent among them. Then, it fol- 
lows that for any E < 0, 

-X 2 a N+ i- n {E) +ui < K n (E) < -X 2 a N+1 - n {E) +u N , 

(30) 

and UN+i-n{E) ^ for n = 1, . . . , N- ln( i, while 

UJl < K n {E) < UJ N , (31) 

for n — A; n d + 1, . . . , A. Therefore, only the first Aj n( j 
eigenvalues of K(E) are ensured to be negative as \X\ goes 
to infinity without regard to the location of {w n }„ =1 . 

Proof : Taking —X 2 S(E) as the unperturbed part of 
K{E), we obtain Eq. for all n. Note that if only 

Ajnd form factors are linearly independent, it holds that 
a m (E) = for m = 1, . . . , N — Aj n d and otherwise does 
not vanish. Then, the assertion is proved straightfor- 
wardly. D 
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FIG. 1: ujz — « n (0) for n = 1, 2, 3 (three solid lines) for 
a three-level system with form factors 13211 . plotted against 
A, and uz — ui, uiz — uii (two dashed lines), and uz (dot- 
dashed line), for reference, where Uz — u\ > ^3 > Uz — W2- 
Three different regions are distinguished, corresponding to the 
number of solid lines satisfying uz — «n(0) > Uz, that is, just 
the number of negative eigenenergies of H, by Theorem lIII.4l 



To illustrate the emergence of the negative eigenener- 
gies, described in Theorem llll . 41 and Proposition lIII.5| let 
us consider the three-level system especially in the case 
where U\ < while w 2 > and u 3 > 0. We also choose 
three form factors, such as 



« n (w)=A 1 / 2 



[1 + (u / K) 2 Y+ n 



(32) 



where A is the cut off constant, and a n is a parameter. 
The form factors described by such algebraic functions 
are often found in various systems involving the process 
of the spontaneous emission of photons from the hydro- 
gen atom [22l l^ | .the photodetachment of electrons from 
negative ions |3, EL EJ ) an d quantum dots [2]|. In the 
calculation depicted in Figs. Q and |2 we have chosen 
a set of parameters wi/A = —0.01, L02 / A = 0.01, and 
w 3 /A = 0.02, and ai = 0.0, a 2 = 2.0, a 3 = 1.0. These 
choices for a n guarantee linear independency among v^s 
, so that A^ind = 3. 

Figure n shows w 3 — k„(0) for n — 1,2,3, changing 
A from 0.1 to 10.0, and ^3 — u\, lu 3 — cj 2 j (two dashed 
lines) and L03 (dot-dashed line) for reference. The latter 
satisfy the relation that 073 — 071 > 073 > 073 — 07 2 > 0. 
One may recognize three different regions in this figure: 
for small A < 0.2, one inequality 073 — fti(0) > 073, i.e., 
ki(0) < 0, holds. In the next region 0.2 < A < 1.0, 
two inequalities, 073 — ki(0) > 073 and 073 — k 2 (0) > 073, 
hold. For A > 1.0 the last region, three inequalities, 
^3 — n n (0) > 0J3 for all n = 1,2,3, are satisfied. There- 
fore, according to Theorem llll. 41 one sees that one, two, 
and three negative eigenenergies of H exist in the first, 
second, and third regions, respectively. It is worth not- 
ing that the appearance of the negative eigenenergy in 
the first region merely occurs from the fact that lu± < 
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FIG. 2: Ldz — K n (E) for n = 1, 2, 3 (three solid lines) for 
a three- level system of the form factors (1321 . uz — E (short- 
dashed line), and 073 — and 073 ~ U2 (two dashed lines). 
We plot them in A = 0.1, 0.7, and 10.0, in (a), (b), and 
(c), respectively, choosing the parameters as 07i/A — —0.01, 
U2/A = 0.01, and Uz/A = 0.02. (a) For a relatively small A, 
only ujz — Ki(E) intersects ujz — E at E/A ~ —0.02, which 
is predicted by Theorem 1111.41 Thus there is one negative 
eigenenergy. One also sees that uiz — «i (E) and u>z — «2 {E) 
still lie closely above the asymptotes u>3 — u\ and Uz — U2, 
respectively, (b) Both luz — ni(E) and 073 — ^{E) intersect 
uiz — E in the vicinity of E = —0.3 and 0.0, respectively, so 
that there are two negative eigenenergies. (c) All ^3 — K n {E) 
for n = 1, 2, and 3 intersect uz — E, and thus three negative 
eigenenergies exist. In this figure, only two intersections for 
n — 2 and 3 are depicted. 
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(see, the latter part of Theorem IIII.4|) . whereas that in 
other regions could be understood as a strong-coupling 



effect (Proposition IIII . 

Figure[2shows three curves of u!^ — K n (E) for n = 1,2,3 
(three solid lines) and us^ — E (short dashed line) , plotted 
against E. An intersection of the former and the latter 
means an emergence of a negative eigenenergy. We also 
plot the asymptotes a->3 — lo\ and L03 — 0J2 (two dashed 
lines), to which 0^3 — Ki(E) and 0J3 — K2(E) are close 
from above as E — > —00, respectively (see, Lemma llIL3ll . 
Figures El (a), 12 (b), and |2 (c), are in the cases where 
A = 0.1, which belongs to the first region, A = 0.7, of the 
second one, and A = 10.0, of the last one, respectively. 
See Fig. Q It is seen in Fig. |2 (a) that W3 — E intersects 
Ki(E) only, so that there is one negative eigenenergy. 
In Fig. (b), one distinguishes the two intersections 
between U3 — E and U3 — n\ (E) , and between 103 — E and 
w 3 — K 2(E). Thus two negative eigenenergies appear. 
The intersection between the latter pair still lies around 
E = 0.0. In Fig. [5] (c), where a relatively large A was 
chosen, lo^ — E finally intersects all three lines, — n n (E) 
for n — 1,2, 3, which tells us three negative eigenenergies 
exist. 



IV. ABSENCE OF BOUND-STATE 
EIGENENERGY INSIDE THE CONTINUUM 

Let us next examine the nonnegative-eigenvalue prob- 
lem for Eqs. © and J7J. In this case, the normalization 
condition © does not hold automatically, unlike the case 
where E < 0, because of a possible divergence of f(uj) at 
u> = E. Before going to the iV-level case, let us first 
observe the single-level one. Except in the trivial case 
where c\ = 0, the condition @ for an eigenvalue E > 0, 
if any, imposes the nontrivial condition or constraint that 



Vi(E) = 0, 



(33) 



where we assume some extent of the smoothness of v\{uj) 
|26|. Then, f(ui) = —\ciVi(uj)/(uj — E) is ensured to be 
square integrable, and Eq. JfJJ) reads 



A^ 



E 



■dcu = E. 



(34) 



To find the solution E of Eq. (|34|l . one may attempt 
to interpret it as an intersection between the left-hand 
and the right-hand sides, as in Eq. (|12J) . However, this 
approach seems impossible at first, because the left-hand 
side of Eq. (|34|) is not well defined for a general E ex- 
cept such points satisfying Eq. (|33() . This matter can be 
solved by alternatively considering the following equa- 
tion, 



ui - X 2 P 



E 



•dcu = E, 



(35) 



that is obtained from Eq. (|34|) by replacing 

with its principal value P J °° dcu. In this case, 

the left-hand side can make sense for a general E, and 
we can treat E as an independent variable. If we find 
the solution E of Eq. (|35|) . and furthermore if it satis- 
fies Eq. I|33(l . it becomes a true solution of the original 
equation (|34|l . Indeed, in such a situation, we have that 
joo \j^f dLU = p l^fdu, and thus Eq. G3 just 
reproduces Eq. (|3"4*|) . 

In the iV-level cases, the condition © for an eigen- 
value E > 0, if any, can be translated into the equivalent 
condition for both the coefficients {c n }^ =1 and E, that 



J2cnVn{E) =0. 



(36) 



Under this condition, we can safely substitute Eq. (JSJ) 
into Eq. (jSJ. However, similarly to Eq. I|35|l . we consider 
the alternative equation in the iV-level cases as 



JV 



[^n<W - A 2 ZW(£0]c„' = Ec„ 



(37) 



n' = l 

for n = 1, 2, . . . , N, where 



D nA E) := P r d., (38) 

Jo lo — E 



which are the components of the hermitian matrix D{E) 
defined for all E > 0. One sees that Eq. ll3*7|l has 
the same form as Eq. <(T6T) . except the point that S(E) 
(E < 0) is replaced by D(E) (E>0). Then, we can im- 
plement a formulation in the matrix form, just as in the 
preceding section. In fact, the solutions of Eq. I|37|l can 
be connected with those of Eq. © under the condition 
p)|> . We first note that 



P 



djj 



N 

= E / 



uj-E 



uj-E 



dw, 



(39) 



which is always valid for all E. Then, substituting this 
relation into Eq. I|37|) . we have 



w„c„ - X 2 P 



LO- E 



dcu = Ec n , (40) 



for n = 1, 2, . . . , N. For a comparison, see Eq. again. 
Therefore, if the solutions E and {c„}^ =1 of Eq. ij^Hjl . 
i.e., Eq. lpT7|). satisfy the condition JSSJl, Eq. ijtDJl can 
reproduce Eq. (|HJ), so that the solutions of Eq. I|40() 
become the true ones of Eq. ©. 

Our procedure for finding the coefficients {c n }^ =1 and 
the nonnegative eigenvalue E of H that satisfy Eqs. © 
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and J7J consists of the two steps: we first solve Eq. (|37|l . 
and then we check whether the solutions satisfy the con- 
dition (35} • For a later convenience, we introduce the 
hermitian matrix K(E) for E > whose components are 
defined by 

K nn , (E) := u n 6 nn , - X 2 D nn , (E). (41) 

Then, the existence of a nontrivial solution of Eq. 1|37|) 
is ensured if and only if there exists a nonnegative E to 
satisfy 

K n (E,X)=E, (42) 

for a certain integer n, where {n n (E, A)}^ =1 are the 
eigenvalues of K(E), arranged in increasing order. To 
summarize again, if an eigenvalue of K(E) is E, then it 
is an eigenvalue of H, provided that it also satisfies the 
condition Q36|). 

It is worth noting that the condition seems not 
necessarily to require the existence of a zero of v n (ui), 
unlike the single-level case of Eq. I|33(l . However, the fol- 
lowing statement means that if v n (u> n ) ^ for all tu n > 0, 
the weak-coupling condition results in no positive eigen- 
value of H strictly. 

Theorem IV. 1 Suppose that Hq has N + positive 
eigenvalues without any degeneracy, and each v n (u>) 
is an L? -function of the form, v n (uj) — uj Pn f n (uj), 
where p n > and f n {^) is a C 1 -function in 
[0,oo). Furthermore, it is assumed that there is some 
Sq > such that sup w>5o \v^(u})v n ' (lu)\ < oo and 
sw Pu>s \d[ v n( LU ) v n'( UJ )]/dw\ < oo for alln andn' . Then, 
if X is sufficiently small but not zero and the condition 
that v n {oj n ) ^ for all n > N — N + + 1 is satisfied, H 
has no positive eigenvalues. 

Proof : Under the assumption that E > 0, we first 
consider the eigenvalue problem 

N 

{E,X)c m , (43) 

i=i 

for i = 1, 2, . . . , N, where {c n i\f =1 is the normalized 
eigenvector corresponding the n-th eigenvalue k„(E,X) 
of K{E). Then, by Theorem 4.3.1 in Ref. [U, one sees 
that 

\K n (E,X)-w n \ < X 2 max{\d 1 (E)\,\5 N {E)\} 

= A 2 ||Z?(_E)| < A 2 sup ||D(.E)||, (44) 

for all n, where 5 n (E) is the n-th eigenvalue of 
D(E). Note that from the assumption in the theo- 
rem and the propositions in the appendix, it holds that 
sup B>0 ||D(.E)|| < oo. Therefore, by choosing A so that 
|A| < A a , we have the fact that \n n (E, A) — uj n \ < R a 
for all E > and all n, and in particular K n (E,X) > 



for all E > for all n > N - N+ + 1, where A Q = 
{Ra/sup^WDiE)]]) 1 / 2 and 

R a = min <^ Wjv-w + +i/3, min{|w„ - w m |/3 \n^m} \ . 

I n,m I 

(45) 

The latter means that k 71 (E,X) only for n > N — 
N + + 1 becomes a candidate for positive eigenvalue 
of H . Note that km~n + cannot be such a candidate 
even if lon-n + — 0. Because in such a case, putting 
X b = {R b /sup E>0 \\D(E)\\)^ 2 , wc find from Eq. O 
that for |A| < X b , \k N - N+ (E,X)\ < R b for all E > 0. We 
here choose such a R b as to satisfy that D(E') > for all 
positive E' < R b . Existence of such an R b is ensured by 
Eq. 1A.2[I in Proposition^ Then, from Theorem 4.3.1 
m Ref. 20] again, we have the estimation that 

-X 2 6 N (E) < k n ^ n+ (E,X) < -X 2 5 l {E) < 0, (46) 

for all E < R b . Hence, we conclude that if |A| < X b , it 
holds that k N - N+ {E, A) < E for all E > [iffi . 

However, we can show that if we choose A sufficiently 
small, any such a K n (E, A) and eigenvector J2i c ni\i) can- 
not satisfy Eq. I|3tj|) . no matter how well they satisfy Eqs. 
(33 and O- To this end, let us look at Eq. (35J, which 
is rewritten as 



N 


2 


2V 2 


i=l 




P n {E, A) £ V*(K n )\i) 



= KWI 2 



N N 

+ ^2(i\v i ( Kn )[P n (E, A) - \n)(n\] ]T 10, (48) 

i=l i' = l 

where P n (E,X) denotes the projection operator associ- 
ated with the n-th eigenvalue K n (E,X). One sees that 
the first term on the right-hand side of Eq. I|48|) behaves 
as 

lim|<( K „(i?,A))| 2 = KK)| 2 , (49) 

for all E > uniformly, because of Eq. 144fl . From the 
assumption of the theorem, \v* t (w„)| 2 does not vanish. 
For the second term on the right-hand side of Eq. i|48[) . 
we can use the result on the perturbation of the projec- 
tion operator [2jJ, which leads to the fact that 

oo 

P n (E,X) = \n)(n\+ £A 2 'PW, (50) 
i=i 

with 

(Ko-cr'iDmKv-cr'Ydc, (51) 

where r„ is the closed positively-oriented circle around 
£ = u> n with radius min m (^„){|w„ — uj m \/3}. Series <|50[) 
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is ensured to converge uniformly for all A such that |A| < 
min{A a , A;,}, because 



sup |A| 2 p(£;)||||(Ko -Cr 1 !! < A 2 /A 2 <l, (52) 
Cer„ 



where 

A n = 



min {\uj„ - w TO |/3} / sup ||-D(-E)|| 

m{=jtn) I E>0 



1/2 



(53) 



From the assumption of no degeneracy among {oJ n }n=i 
and the discussion after Eq. (|44l) . for such a A, all r n 's 
are disconnected from each other, and there should be 
only one eigenvalue of K in each circle. This leads to 
dim[P n (E, \)C N ] = dim[\n)(n\C N ] = 1, so that A = 
is not an exceptional point 29]. It is worth noting that 
A„ does not depend on E. Thus, the second term on the 
right-hand side of Eq. (|48|l is estimated as 



N 



N 



J2mvM[P n (E,X)-\n)(n\}J2^( K nm 



< 



< 



N 



i=l 
N 



Y su p KM I 



\\P n (E,X)-\n)(n\\\ (54) 
(A/A„) 2 



1 - (A/A„)2 



0,(55) 



as A — > 0, for all E > uniformly, where it was used 
that sup Cer?i $ r \\(K - C) 1 1 1 1 <^CI < 2tt. Eq. gHJ with 
the results (|49|l and l|55|l means that Eq. I|36[l is never 
satisfied for sufficiently small A with |A| < min{A a ,A{,}, 
even if K n (E, A) = E holds. D 

It is worth considering the opposite condition that 
ii n (u„) = 0. In this case, we could infer the existence of 
an eigenvalue inside the continuum, from the decay pro- 
cess arising from the pole z p , n - Indeed, if we recall the 
explicit form of the decay rate [3(j , if the opposite con- 
dition holds, the decay rate comes small so that a much 
slower decay occurs. Then, one may associate such a be- 
havior with the presence of a bound state [3l| , though 
it is not obvious whether this pole actually becomes an 
eigenenergy of H. 

Let us now evaluate an explicit value of A for which 
there is no positive eigenvalue of H. Under the as- 
sumption of the analyticity of v n , one sees that if |A| < 
min{A a ,Afc}, Eq. (|49|l is rewritten by using Eqs. 144|) 
and lf55)) as 



\\v*(K n (E,X))\ 2 - |<(w„)| 2 
d\v n (uj)\ 2 



< sup 

\u)—Lj n I <R a 



A 2 

= tt min {|w„ 



w m |/3} sup 

LU>0 



\K n (E,\)-u n \ (56) 
d\v n {u)\ 2 



(57) 



Therefore, by setting Eqs. (JSSJ) and (JEJ hito Eq. (|4Tjl . 
the left-hand side of Eq. 147|) is ensured to be positive, 



and no positive eigenenergy of H exists, providing that 
A is chosen to satisfy the N + + 1 inequalities, 



|A| < min{A Q , A & }, 



(58) 



and 



A 



K(w„)| > -y min {\w n - w m |/3}sup 



A 2 m(=£n 



N 



Y sup \Vi(u)\ 



UJ>0 

2 



d\v n (uj)\' 



A 2 /A 2 



1 - A 2 /A; 



-, (59) 



for n = N - N + + 1, . . . , N. By solving Eq. for A 
explicitly, Eqs. l(55|> and are reduced into the single 
inequality 



|A| < min{A a , A h ,A 



N-N++1) ■ ■ ■ , 



Ajv}, (60) 



with 




Un + fin + In ~ V (<*„ + (3 n + 7„) 2 - 4a„/?„] 

(61) 



where a n = |<(w„)| 2 , (3 n 
w m |/3} sup a)>0 \d\v n (uj)\ 2 / duj\, 



and 7„ 



i=l SU PI 



U~U> n \<R c 



In order to demonstrate Theorem lIV.ll we apply it to 
the spontaneous emission process for the h ydr ogen atom 
interacting with the electromagnetic field [23. We sup- 
pose that |n) is the product state between the (n + In- 
state of the atom and the vacuum state of the field, and 
also \u>) the product state between the ls-state of the 
atom and the one-photon state. Then, an initially ex- 
cited atom is expected to make a transition to the ground 
state by emitting a photon. We treat the atom as a four- 
level system composed of the ground state and the three 
excited states: the 2p, 3p, and 4p state. The form fac- 
tors corresponding to the 2p — Is, 3d — Is, and Ap — Is 
transitions were obtained as follows [13, 123, 1231 , 



iA 



1/2 (lo/A-l) 1 / 2 
1 [l + Cw/AO^' 

1/2 (aVA 2 )V 2 [l + 2(a,/A 2 ) 2 ] 
1 128V2[1 + (^/A 2 ) 2 ]3 ' 

i54V3A} /2 (w/A 3 ) 1/2 

45 + 146(w/A 3 ) 2 + 125(w/A 3 ) 4 
X 15625[1 + (w/A 3 ) 2 ] 4 



(62) 
(63) 

, (64) 



where Ai = 8.498 x 10 18 s^ 1 , A 2 = (8/9)Ai s _1 , and 
A 3 = (10/12)Ai s _1 are the cut off constants. One sees 
that these form factors satisfy all conditions required in 
Theorem IIV.1I The coupling constant is also given by 
A 2 = 6.435 x 10 -9 . The eigenvalues of H are given by 
Lo n = ffi[l - (n + 1)~ 2 ] with Q = 1.55 x lO^s" 1 , all 
of which are embedded in the energy continuum. The 
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Hamiltonian £[J is then derived under the four-level ap- 
proximation (i.e., N = N + = 3) and the rotating- wave 
approximation. The various parameters are numerically 
obtained as follows: R a = \oj 2 — W3I/3 = (7/324)0, 
su PjB>0 \\D(E)\\ = -61(E) = 11.332Ai at E = 0.6145Ai, 
A? = 5.45 x 10- 8 fi/Ai, X\ = A§ = X 2 a = 1.91 x l(r 3 fi/Ai, 
ai = 1.82 x 10" 3 Ai, a 2 = 4.87 x 10- 4 Ai, a 3 = 1.99 x 
10" 4 Ai, Pi = 6.17 x 10- 2 fl, p 2 = 4.87 x 10- 3 n, p 3 = 
1.88 x 10- 3 Q, 7l = 2.45 x 10" 3 Ai, 72 = 3.04 x 10" 3 Ai, 
73 = 2.45 x 10~ 3 Ai, from which Eq._ JJHJ reads X\ = 
4.18 x 10~ 6 , A 2 , = 5.01 x 10~ 7 , and X 2 = 2.14 x 10~ 7 . 
Then, it follows that 



mm{X 2 a ,Xl\m} = A§ > A 2 , 



(65) 



and thus Eq. (|60|l holds. This conclusion indicates that 
the intrinsic values of the parameters characterizing the 
system does not allow any bound state. In fact, we have 
not observed any such state. It is worth noticing that the 
upper bound estimated in Eq. (|65|l is dominated by the 
factor A3, roughly speaking, the minimum level-spacing 
over the maximum cut off constant. 



V. CONCLUDING REMARKS 

We have considered the eigenvalue problem for un- 
stable multilevel systems, on the basis of the JV-level 
Friedrichs model, where the eigenenergies are supposed 
outside or possibly inside the continuum. The outside 
case is essentially determined by the location of the dis- 
crete level u> n of the free Hamiltonian and the strength of 
the coupling constant A. If ui n lies outside the continuum, 
the corresponding eigenvalue always lies below uj n . If ui n 
lies inside the continuum, by choosing a A large enough 
the eigenvalue originating from w„ can emerge from the 
continuum. Such behaviors are similar to those seen in 
single-level cases, however, this is not the case if the form 
factors v n are linearly dependent. On the other hand, we 
have shown the absence of the eigenvalue lying inside the 
continuum in the weak coupling cases, under the condi- 
tion that v n (uj n ) ^ if u n lies inside the continuum. This 
statement is just an extension of Lemma 2.1 in Ref. [l4| . 
where only identical form factors were considered, and 
the upper bound for |A| required in the lemma was not 
estimated. We have evaluated this upper bound in our 
case, which proves to be proportional to the minimum 
level-spacing over the maximum cut off constant. Hence, 
comparing this value with the actual A, one can check 
at least the absence of the eigenvalue, even in the case 
that one cannot evaluate the reduced resolvent explicitly. 
At first sight, the normalization condition, i.e., Eq. (|36|l . 
seems not necessarily to require the zeros of the form fac- 
tors for a presence of a bound-state eigenenergy inside 
the continuum, though it is misplaced in weak-coupling 
regimes. However, we still do not have a definite an- 
swer to this matter in other coupling regimes where the 
multilevel effect may allow a presence of a bound-state 



eigenenergy inside the continuum without zeros of the 
form factors. 
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APPENDIX 

In this section, we present the two Propositions ["71 and 
f2l The former and the latter state that the behavior of 
the energy shift D(E) at small and large energies is quite 
regular without any divergence, respectively, under some 
form-factor conditions that are often satisfied by actual 
systems. 

Proposition. 1 Suppose that the function r)(ui) belong- 
ing to L 1 ([0, 00)) is of the form 



r](u) := Lu p r(uj), 



(A.l) 



where p > and r(u) is a C 1 -function defined in [0, 00). 
It then holds that r)(u))/uj S L 1 ([0,oo)) and 



— — -dw = lim f — — ^—diu = limP f — — — \-dw. 



-Eto J uj - E Eio J Q uj — E 

(A.2) 



Proof : From the proof of Proposition 3.2.2 in Ref. 
|l6j | , the principal value of the integral on the right-hand 
side is written by the absolutely integrable function as 
follows 



P 



auj = 

W-E ./o 



T](bj) - r}(E)(f S (oj - E) 
oj-E 



doj, 
(A.3) 

for all E > 0, where (ps(uj) is a C^-function with support 
[—5,5] (0 < V5 < E), even with respect to the origin, 
and such that <ps(0) = 1. In the following, we choose 
ip s {uj) = exp[l - 1/(1 - (uj/S) 2 )} for lo £ (-6,6) or 
otherwise, and 6 — E/2. On the other hand, since from 
the assumption (jA.lll T ](u)) [to is absolutely integrable, the 
first equality in Eq. 1A.2I) is obvious. Therefore, it is 
sufficient to show that 



lim 

£10 



T}(w) rj(ui) - r/(E)ip s (uj - E) 



E 



0. 

(A.4) 
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Note that the above integrand can be rewritten as 
r)(u)) r](u) - T)(E)<p s (uj - E) 



= -E 



T)(uj) 



uj-E 

r)(E)<p 5 (u - E) 



+ 



uo(uj — E) 
7]{E) m {u-E) 



E 



uj-E 

■qjuj) - T](E)ips(uj - E) 
w(w - E) 



(A.5) 
(A.6) 



Let us first consider the case where ui € / := (0, E/2] U 
[3E/2,oo). Then, since (ps(uj — E) = 0, we can use Eq. 
(|A.5|I to estimate the integrand: 



E- 



T){uj) 



uj(uj - E) 



< 2 



(A.7) 



where the right-hand side is absolutely integrable and 
independent of E. Furthermore, it follows that 
liniEjo Exi(oj)rj(uj)/[uj(uJ — E)] = for every ui S (0, oo), 
where x/( w ) — 1 (w £ /) or (w ^ J), being the char- 
acteristic function. Thus, by the dominated convergence 
theorem, we can see that 



lim 

EJ.0 



E/2 



E 



T)(UJ) 



3E/2 



uj{lj - E) 



dw = 0. 



(A. 



For ui £ (E/2,3E/2), we can use Eq_ (O) . The in- 
tegration of the first term of Eq. (|A.6(1 is estimated by 



3E ^ r,(E)M^-E) 

E/2 ^ 

rj(E) r^/2 



dui 



- ~lk7K I ips(u - E)duj = r)(E) ip 1 (x)dx^>Q, 

^1 1 JE/2 J-l 

(A.9) 

as E [ 0. The second term of Eq. (|A.6|I is also estimated 
by 

\v(u) - rj{E)¥>s(u - E)\ 

< \ V (uj) - V (E)\ + \ V (E)\\1 - ip s (u> - B)|.(A.10) 

The integration of the first term on the right-hand side 
right-hand side of the above is evaluated as 



E/2 



LO\LO - E 



<(ln3)£ sup |t/(w)I 

E /2<u<3E /2 



(A.ll) 



< (ln3)£ 



p J BP- 1 max{(i)P- 1 ,(|)P- 1 } sup |r(w) 



( ¥P) sup 1/(07)1 

' / uj£[0,3E/2] 



wG[0,3E/2] 

as E | 0, (A.12) 



where the prime on rf (ui) implies the differentiation of 
7](u>) and so on. The integral corresponding to the last 



term on the right-hand side of Eq. (|A.10|) is also esti- 
mated as 

3E /* e \ v (e)\\i-M"-e)\ ^. 

E/2 U)\U)-E\ 

< (hi 3)E\ V (E)\ sup \<ff s (u-E)\ (A.13) 

E/2<ui<3E/2 

= 2(ln3)\ V (E)\ sup |^(x)| (E j 0). (A.14) 

kl<i 



Thus, we can obtain 

lim r\^-^i { r E) ^ = 0- (A-15) 
EioJ E/2 uj(uj-E) 



Equations (|A. 8|) , (| A. 9|> , and (|A.15fl mean the completion 
of the proof of HA.4J1 . □ 

Proposition. 2 Suppose that the function rj(uj) be- 
longs to L 1 ([0, oo)) H C ([0, oo)), and satisfies that 
su Pw>,5 D \v( u )\ < 00 and s ^Pu>8 W{u)\ < oo for some 
Sq > 0. Then, 



sup 

E>S 



P 



uj-E" 



du 



< oo. 



(A.16) 



Proof : To examine this integral, we use the expression 
(|A.3(I and divide the interval [0, oo) into Is.e = [E — 
S.E + 8] and Is,e — [0, oo) \ Is t E, again, where we assume 
So > S > 0. In the latter interval, it is estimated that 
Xi ^(ujM^/(uj-E)\ < \tj(uj)\/S e^ao.oo)). Then, 



sup 

E>S 



y 7 — (ui) duj 



1 f°° 

< -r / \r)(uj)\duj < oo. 
6 Jo 



In the former interval, the integrand in Eq. (| A. 3|> is eval- 
uated as 



77(0;) - t)(E)lps(uj - E) 



uj-E 



which results in 



sup 

-E><5 



< 2(5 



< sup \rf(w)\ 

+ W(E)\ sup \<pe'((A\m 

|o>|<<5 



r)(uj) - r](E)ips(uj - E) 



uj-E 



sup \r,'(E)\ + sup \r)(E)\ sup \<p s \w)\ 

E>S E>S a \<»\<8 



< CO, 

(A.19) 



where we used the assumption for rj(u>) in the statement. 
Incorporating Eq. I|A.17|) with Eq. i|A.19|) . Eq. I|A.16|) is 
obtained. D 
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